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Abstract 


This  report  summarizes  same  of  the  research  carried  out  under  contract 

F  61052-69-C-0020  during  the  year  1st  September  1968  to  31st  August  1969. 

laminar  and  transitional  heat  transfer  rates  were  measured  in  various 

base  type  separated  flows  in  the  shock  tube.  It  was  found  that  maximum  and 

l/2 

average  heat  transfer  rates  are  correlated  by  the  parameter  (hRe^  /l)  for 

models  with  an  initial  boundary  layer  and  by  Re.  for  models  with  leading  edge 

L 

separation . 

The  application  of  the  hot  wire  technique  f^r  measurements  of  the  flow 
field  in  the  supersonic  near  wake  is  studied.  It  is  shown  that  using  the  mea¬ 
surements  of  a  hot  wire  probe  and  a  total  pressure  probe  enable  the  determina 
tion  of  all  the  required  physical  parameters  of  the  dissipative  flow  field  in 
the  near  wake.  The  flow  field  behind  a  blunt  two  dimensional  base  model  is 
measured  at  Kach  number  of  2.25  and  Re^  =  1.5  x  10^.  Mach  number  and  total 
temperature  profiles  are  determined  and  the  various  zones  in  the  flow  field  are 


identified. 
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I.  INTRODUCTION 

The  present  report  summarizes  some  of  the  research  work  carried  out  under 
contract  F  61052  69  C  0020  during  the  year  1st  September  1968  to  31st  August 
1969.  This  year’s  program  included  an  extension  of  the  heat  transfer  studies 
of  base  type  flows  in  the  shock  tube  which  are  discussed  in  Part  II,  as  well 
as  a  program  of  near  wake  flow  field  study  using  the  hot-wire  anemometer  dis¬ 
cussed  in  Part  III. 

II.  Heat  transfer  rate  measurements  in  base  type  flows  on  the  shock  tube 
Flow  separation  at  supersonic  and  hypersonic  speeds  is  known  to  have 

strong  effects  on  the  heat  transfer  rates  to  body  surfaces  particularly  in 
the  reattachment  regions.  In  general,  it  was  found  that  heat  transfer  rates 
are  decreased  in  the  ’’dead  water"  zone  and  increased  in  the  reattachment  zone 
in  relation  to  flat  plate  heat  transfer  rates  under  similar  flow  conditions. 

The  net  effect  of  this  variation  of  heat  transfer  on  the  heat  transfer  rate  in 
the  separated  flow  region  as  a  whole  is  not  easily  determined  since  in  many 
cases  the  increase  in  the  heat  transfer  rate  in  the  reattachment  zone  is  much 
greater  than  the  corresponding  reduction  found  in  the  "dead  water"  zone.  For 
the  past  several  years,  heat  transfer  rates  in  separated  flows  have  been  stu¬ 
died  at  the  Aerodynamics  Laboratory  of  the  Technion'e  Aeronautical  Department 
supported  in  part  by  SOAR  contracts.  These  studies  were  conducted  in  the 
3"  x  3"  shock  tube  and  have  included  the  following  model  configurations  : 

(1)  Two-Dimensional  Backward  Pacing  Step  (Ref.  l). 

(2)  Axisymmetric  Backward  Pacing  Step  (Ref.  2). 

(3)  Two-Dimensional  Blunt  Base  (Ref.  3). 


(4)  Sharp  Protruding  Leading  Edge  (Ref.  4). 

(5)  Two-Dimensional  Leading  Edge  Separation  Bubble  (Ref.  5). 

In  this  report  a  summary  of  some  of  the  more  significant  results  of  the 
previous  studies  is  presented  and  is  also  correlated  to  determine  the  main 
parameters  affecting  the  data. 

The  shock  tube  has  been  found  to  be  a  very  useful  facility  for  heat  trans¬ 
fer  studies.  Using  thin  (or  thick)  film  pluf  resistance  then:  **  '■■tors  local 
heat  transfer  rates  can  be  measured  to  an  accuracy  of  better  than±l^  and  even 
up  to  ±  5$  •  These  are  as  good  as,  nr.d  in  many  cases,  even  better  accuracies 
than  can  be  obtained  in  heat  transfer  measurements  in  conventional  wind  tunnels. 
The  feasibility  of  measurements  of  heat  transfer  in  separated  flows  in  the  shock 
tube  and  the  shock  tunnel  was  discussed  by  a  number  of  investigators  (Ref.  6,  7, 

and  8).  Thus,  based  on  measurements  in  the  shock  tube  our  previous  results 

/,v  / 

Indicated  that  laminar  heat  transfer  rates  fre  separated  flows  (with  an  attached 
boundary  layer  ahead  of  separation)  can  be  correlated  using  the  parameter 

hRe. V  /L.  This  parameter  is  proportional  to  the  ratio  of  step  height  to  boun- 

L 

dary  layer  thickness  at  separation,  ,  and  was  found  to  correlate  the  laminar 
separated  region  pressure  field  as  well.  It  is  interesting  to  note  that  some 
measurements  of  the  heat  transfer  rate  distribution  on  various  back-ward  facing 
stsp  models  in  a  continuous  supersonic  wind  tunnel  reported  in  Ref.  9  support 
this  correlation.  Comparison  of  the  shock  tube  data  and  the  wind  tunnel  measure¬ 
ments  reported  in  Ref.  9  are  included  in  this  paper.  Although  the  flow  condi¬ 
tions  in  the  two  experiments  are  very  different,  particularly  the  large  diffe¬ 
rences  in  stagnation  to  wall  temperature  and  the  state  of  the  boundary  layer  of 
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these  experiments,  the  heat  transfer  rates  in  both  exp  rrc.otu  ur  •  fo  *,  1  u,  bn 

l/2 

related  by  the  parameter  hfie.  /L. 

The  experimental  apparatus  and  results  are  discussed  in  .  m  in  iiof.  10. 
The  main  results  are  discussed  herein.  The  variation  of  the  loca.  heat  trans¬ 
fer  rate  along  the  model  length  as  shown  in  Figs.  3*  4,  5  and  6  for  va:  ciu, 

Mach  numbers  indicates  that,  for  a  given  model,  there  is  only  slight  van one 
in  the  position  of  the  beginning  of  the  reattachment  zone  and  the  point  of  tin 
maximum  local  heat  transfer  rate  as  shock  Mach  varies  from  5  to  10.  Detection 
of  any  movement  of  these  positions,  however,  was  limited  by  gage  size  and  loca¬ 
tion.  Of  importance  is  that  all  the  curvea  exhibit  qualitatively  aimilar  va¬ 
riations  of  local  heat  transfer  rates.  This  similarity  of  the  heat  transfer 
rate  distributions  is  found  in  spite  of  large  variations  (flow  Mach  numbers 
and  stagnation  to  wall  enthalpy  ratios)  as  well  &b  model  geometry.  With  the 
exception  of  the  leading  edge  separation  bubble  on  the  blunt  nosed  body,  the 
heat  transfer  distribution  over  the  sharp  leading  edge  model,  the  backward  fac¬ 
ing  step  and  the  axisymaetric  step  models  definitely  show  a  zone  of  low  heat 
transfer  rate  in  the  "dead  water"  zone.  Following  this  is  the  reattachment 
zone  as  signif  , ed  by  a  sharp  rise  ir.  the  heat  transfer  rate.  It  is  expected 
that  following  reattachmert  the  heat  transfer  rates  would  approach  the  flat 
plate  values.  In  fact,  many  of  the  curves  exhibit  this  tendency,  however,  in. 
some  cases,  the  gages  available  on  the  model  were  not  sufficient  to  show  this 
trend. 

The  extent  of  the  low  heat  transfer  rate  zone  in  ‘.he  "dead  wa U  r"  :\- 
is  fcur.i  to  be  about  1  to  1.5  step  heights  from  the  two-dimensional  and  «xi- 
symmetric  backward  facing  steps  (Figs.  3  and  4)  while  it  ©xtor.du  to  a*  *  ■' 
step  heights  for  the  sharp  lending  edge  model  (Fig.  r).  Tr.e  mix'..-:  -  •' 
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transfer  rate  in  the  reattachment  zone  is  found  to  occur  at  about  4  to  6 
step  heights  for  the  two-dimensional  and  azisymmetric  backward  facing  steps 
(Figures  3  and  4)  and  at  about  8  stey  heights  for  the  sharp  leading  edge 
model  (Fig.  5).  For  the  leading  edge  bubble  case,  closure  of  the  bubble  and 
reattachment  is  indicated  at  about  3  nose  heights  behind  the  leading  edge 
(Fig.  6). 

The  results  of  wind  tunnel  measurements  of  the  heat  transfer  rate  distri¬ 
bution  behind  a  backward  facing  step  reported  in  Ref.  9  are  also  included.  It 
is  seen  that  the  wind  tunnel  data  seems  to  be  extended  so  that  the  peak  values 
of  heat  transfer  rates  are  obtained  after  about  1<;  to  15  step  heights  while 
the  shock  tube  data  indicates  that  this  occurs  after  about  4  to  6  step  heights. 
On  the  other  hand  the  results  of  Ref.  11,  which  are  also  obtained  in  the  wind 
tunnel,  show  peak  values  at  about  2  to  3  step  heights.  This  does  raise  the 
question  of  whether  or  not  these  differences  are  dependent  on  errors  due  to  the 
various  experimental  methods  used.  The  heat  transfer  measurements  in  the  wind 
tunnel  require  measurements  of  a  number  of  physical  parameters  which  must  be 
corrected  and  are  sensitive  to  errors,  while  the  measurement  of  heat  transfer 
in  the  shock  tube  involves  a  direct  measurement  of  a  single  precalibrated 
transducer  .  Most  separated  flow  measurements  do  indicate  reattachment  at  2  to 
4  step  heights  therefore  provide  an  additional  support  for  the  shock  tube  mea¬ 
surements  . 

The  ma;  imum  heat  transfer  rate  at  reattachment  is  correlated  by  the  use 
of  the  non-dimensional  parameters  as  follows  : 

(l)  in  cases  of  an  initial  boundary  layer  at  separation 

,  .  Adah/VlO"  Vp_ 


(2)  in  cases  of  zero  initial  boundary  layer  at  separation 


(e.g.  -  sharp  leading  edge); 


The  empirical  values  of  the  parameters  A  and  n  and  B  and  m  for  the 
various  geometries  are  presented  in  Tables  1  and  2,  respectively  and  in  Pig.  7. 
The  maximum  heat  transfer  rate  values  obtained  in  the  wind  tunnel  study  of  Ref. 

9  are  found  to  be  in  surprisingly  good  qualitative  agreement  with  the  shock 
tube  data  in  spite  of  previous  remarks  and  of  the  extreme  differences  in  stag¬ 
nation  to  wall  enthalpy  as  well  as  state  of  flow.  This  seems  to  suggest  that 
the  heat  transfer  in  laminar  separated  flow  is  not  too  sensitive  to  variations 
in  flow  Mach  number  (at  least  in  the  range  of  Kach  number  1.5  to  3.0)  and  stag¬ 
nation  to  wall  enthalpy.  The  dominating  parameter  is  found  to  be  the  value  of 
hRe^/  /L  for  cases  with  initial  laminar  boundary  layer  thickness  or  Re^  for 
cases  of  zero  boundary  layer  at  separation.  additional  important  parameter 
seems  to  be  the  state  of  the  boundary’  and  mixing  layers  and  the  relative  posi¬ 
tion  of  transition.  It  is  suggested  in  our  measurements  and  in  those  of  Baker 
and  Martin  (Ref.  ll)  and  actually  demonstrated  in  the  measurements  of  Sanford 
and  Ginoux  (Ref.  9)  that  when  transition  appears  ahead  of  reattachment  very 
high  values  of  heat  transfer  rates  at  reattachment  are  experienced.  In  the 
completely  laminar  separation  case  the  heat  transfer  rates  are  relatively  low 
and  even  at  reattachment  the  maximum  values  are  not  much  larger  then  the  corres- 
ponding  attached  flat  plate  values.  In  the  completely  turbulent  case  maximum 
values  are  in  order  of  2  to  4  times  the  attached  flat  plate  value  as  found  by 
Naysmith  (Ref.  12)  and  Thomann  (Ref.  13).  However,  in  the  transitional  flow 


cases,  which  we  assume  occured  in  the  shock  tube  tests  at  the  higher  Reynolds 
numbers,  extremely  high  peak  heat  transfer  rates  of  the  order  of  5  to  7  times 
the  laminar  flat  plate  values  were  detected.  The  measurements  in  Ref.  9  tend 
to  support  the  explanation  that  these  high  heat  transfer  rates  are  associated 
with  the  movement  of  the  transition  point  ahead  of  the  reattachment  position 
into  the  mixing  layer.  Therefore  the  present  correlation  may  be  taken  to  cover 
both  the  cases  of  completely  laminar  separation  and  transitional  cases  where 
the  boundary  layer  is  laminar  ahead  of  separation  and  transition  occurs  at  a 
position  between  separation  and  reattachment. 


III.  Application  of  hot  wire  anemometer  for  eupersonic  near  wake  atudiea 

Determination  of  the  physical  parameters  in  the  near  wake  flow  require 
the  measurements  of  at  least  three  independent  local  physical  quantities.  In 
most  previous  wake  flow  measurements  one  physical  parameter,  e.g.  -  the  total 
temperature  has  been  assumed  tc  be  constant  in  the  flow  field.  Measurements 
using  the  hot  wire  probe  and  a  total  pressure  probe  enable  the  determination 
of  all  the  required  physical  parameters  of  the  dissipative  flow  field  in  the 
near  wake.  Due  to  its  relatively  small  dimensions,  fast  response  and  high 
signal  to  noise  output  the  hot  wire  anemometer  is  particularly  suitable  for 
flow  measurements  in  the  wake.  This  part  summarizes  the  method  of  measurements 
and  data  analysis  for  such  measurements  in  the  turbvlent  near  wake  behind  a  two- 
dimensional  blunt  base. 

The  measurements  are  conducted  in  the  10  cm  x  10  cm  supersonic  wind-tunnel 
at  Mach  number  2.25  and  Reynolds  number  based  on  model  chord  of  1.5x10^.  The 
model  is  a  wedge-flat  plate  with  a  blunt  base  and  spans  the  wind-tunnel  test 
section.  The  hot  wire  probe  is  a  0.005  mm  platinum  coated  tungsten  wire  solder¬ 
ed  to  the  two  support  needles.  The  wire  length  is  about  1  mm  (wire  aspect  ratio 
of  about  200).  Two  probe  geometries  (Pig.  8)  were  tested  in  order  to  check 
the  effects  of  these  supports  on  the  directional  wire  response  in  the  recircu¬ 
lating  flow.  Both  probe  supports  gave  essentially  identical  response,  so  that 
this  effect  is  found  to  be  very  small.  The  block  diagram  of  the  measuring 
system  is  shown  in  Fig.  9.  Each  probe  is  calibrated  to  determine  its  resistivi¬ 
ty  coefficient  and  is  pretested  in  a  calibration  wind-tunnel.  A  detailed  des¬ 
cription  of  the  experimental  system  is  presented  in  Ref.  14. 


The  anemometer  response  is  due  to  heat  transfer  from  or  to  the  probe's 


wire  from  the  flowing  fluid.  The  analysis  of  the  data  required  finding  a 
relation  between  the  heat  exchange  coefficient,  indicated  by  the  wire  Ku.mm-I t 
number,  and  the  local  Kach  and  Reynolds  numbers  of  the  flow.  Bnpirical  corre¬ 
lation  of  available  data  was  presented  in  Ref.  15.  For  a  constant  temperature 
anemometer,  used  in  the  present  inveatigation,  the  probe  response  Is 


Nu 

s 


.  R 

r  r 

TFT 

0 


d( V2  /R^) 
d(R  -Raw  ) 


i  >  0 


(1) 


and  the  relation  for  Nu(Re,M)  is  defined  in  Ref.  15  by  relations  of  the  form 


Nu(Re,K)  «  Nu(Re,-)  Q  (Re,H)  (2) 


and  is  shown  in  Fig.  10.  In  a  dissipative  flow  it  is  convenient  to  use  u  moai- 
fled  correlation  of  equation  (2)  writing  the  Mach  number  as  a  function  of 
Nusselt  and  Reynolds  numbers  shown  in  Fig.  11. 

The  adiabatic  temperature  measurement  is  then  corrected  for  the  Knudsen 
and  Mach  numbers  effects  by  the  relation  (Ref.  15) 

n  =  n  (K,Kn)  (?) 

shown  in  Fig.  12. 

Thus,  the  first  stage  in  the  data  a.ialysis  is  calculating  for  each 

measured  hot  wire  output  the  local  Nu  and  Taw  ,  using  the  proper  relations 

m  n 

represented  by  equations  (l),(2)  and  (5).  These  profiles  are  nomalizod  by  the 
free  stream  values  in  order  to  account  for  probes  errors  and  losses.  This  cal¬ 
culation  is  performed  on  the  Elliot  503  computer  using  the  program  shown  in 


block  diagram  in  Fig. 


13. 
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The  local  Mach  and  Reynolds  number  profiles  are  evaluated  from  the  pre¬ 
viously  determined  Nu  and  Taw  by  an  iteration  procedure  shown  in  a  block 

m  m 

diagram  in  Fig  14.  For  this  calculation  the  local  total  pressure  is  obtained 
from  previously  measured  pitot  pressure  profiles  presented  in  Ref.  16. 

The  hot  wire  output  at  various  wire  overheat  values  is  shown  in  Fig.  15. 

From  these  results  the  profiles  of  the  Nusselt  number  and  adiabatic  temperatures 
are  evaluated.  Profiles  of  local  Macu  number  and  total  temperature  are  then 
calculated  using  the  pitot  pressure  data  and  the  iteration  procedure  shown  in 

Fi er.  14.  typical  nrofiles  of  Nu  and  Taw  and  of  Mach  number  and  T  measured  at 

mm  t 

about  one  bas«  height  behind  the  model  are  shown  in  Figures  16  and  17  respectively. 
A  map  of  the  Mach  number  and  of  the  total  temperature  in  the  near  wake  is  shown 
in  Figures  16  and  19  respectively.  A  ohadow  photograph  of  the  corresponding 
flow  field  is  shown  in  Fig.  20.  Using  these  measurements  the  various  flow  re¬ 
gions  can  be  now  specified  as  shown  in  Fig.  21. 

It  is  found  that  there  is  an  excellent  correspondence  between  the  hot  wire 
measurements  and  flow  field  photographs  in  respect  to  the  positions  of  the  lip 
shock,  trailing  shock,  expansion  fan  region,  etc. 

Furthermore,  these  measurements  do  show  that  the  clearly  seer,  bright  white 
line  emanating  from  the  base  corner  does  correspond  to  the  position  of  peaks  in 
the  total  temperature  profiles.  These  peaks  are  found  to  be  near  the  outer  edge 
of  the  mixing  zone.  This  may  indicate  a  layer  of  high  shear  in  this  part  of  the 


mixing  zone 
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In  the  "deed  water"  tone  the  hot  wire  probe  la  aeneitlve  to  the  absolute 
value  of  the  local  velocity  with  no  regard  to  lta  direction.  It  la  therefore 
Interacting  to  note  that  Mach  number  valuta  of  about  0.4  to  0.6  are  meaaurod 
near  the  base,  lbe  Mach  nunber  variation  on  the  wake  center  line  la  shown  in 
Fig.  22.  Near  the  reattachaent  tone  the  preaent  iteration  procedure  failed  to 
converge  ao  that  no  Mach  nunber  values  are  available  there.  This  effect  may  be 
due  to  either  very  high  ahear  values  at  this  zone  (which  cause  erroneous  probe 
output)  or  due  to  large  flow  inclination  which  results  in  considerable  error 
in  the  total  pressure  measurements. 
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